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Abstract. Recently proposed Transverse Enhancement Model of nuclear effects in Charge Current Quasi- 
Elastic neutrino scattering [A. Bodek, H. S. Budd, and M. E. Christy, Eur. Phys. J. C71 (2011) 1726] is 
confronted with the MiniBooNE high statistics experimental data. It is shown that the effective large axial 
mass model leads to better agreement with the data. 

PACS. 13.15.+g Neutrino interactions - 25.30.Pt Neutrino scattering 
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1 Introduction 

Charge current quasi-elastic (CCQE) scattering is the most 
abundant neutrino interaction in oscillation experiments 
like MiniBooNE (MB) or T2K with a flux spectrum peaked 
below 1 GeV. Its full understanding is important for detail 
neutrino oscillation pattern measurements. 

Under an assumption that the impulse approximation 
picture [T] is valid the CCQE reaction both on free and 
bound nucleons is defined as: 
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with v, v, V^,p and n standing for: neutrino, antineutrino, 
charged lepton, proton and neutron respectively. 

A theoretical description of free nucleon target CCQE 
reaction is based on the conserved vector current (CVC) 
and the partially conserved axial current (PCAC) hypothe- 
ses. The only unknown quantity is then the axial form- 
factor Ga(Q 2 ) for which one typically asssumes the dipole 

form Ga(Q 2 ) — Ga(0)(1 + |p-)~ 2 with a free parameter, 

called the axial mass Ma- 

The aim of CCQE cross section measurements is to 
determine the value of Ma and also to validate a nuclear 
physics input used in theoretical cross section computa- 
tions. There is a variety of approaches [2. including the 
Fermi Gas (FG) model implemented in all the neutrino 
Monte Carlo (MC) events generators. 

Measurements of Ma can use an information contained 
in the shape of the distribution of events in the four- 
momentum transfer Q 2 (strictly speaking in the variable 
Qqei see 0) w hich is sensitive enough for precise evalua- 
tions of Ma- The dependence of the total cross-section on 
Ma gives an additional input: if the Ma value is increased 
from 1.03 to 1.35 GeV for E v > 1 GeV the cross-section is 
raised by ~ 30% (for E v < 1 GeV the increase is smaller) . 



Another interesting option to validate models is to com- 
pare to high statistics double differential (2D) cross sec- 
tion data (muon kinetic energy and scattering angle) on 
carbon provided by the MB collaboration [3J. 

In the past, several measurements of Ma were done 
on deuterium for which most of nuclear physics compli- 
cations are absent. Until a few years ago it seemed that 
the results converge to a value ~ 1.03 GeV 1.4] . There is 
an additional argument in favor of a similar value of Ma 
coming from the weak pion- production at low Q 2 . When 
put together they suggest the value Ma = 1.014 GeV 0]. 
On the contrary, all (with an exception of the NOMAD ex- 
periment) more recent measurements of Ma report much 
larger values (for a discussion see: [5]). 

A theoretical mechanism which can explain the Ma 
value discrepancy comes from the many-body nuclear model 
proposed 10 years ago j6] based on the ideas of M. Ericson 
and developed more recently by Martini, Ericson, Chan- 
fray and Marteau (MEChM model) . The model predicts a 
large contribution to the muon inclusive CC cross section 
from elementary 2p-2h and 3p-3h excitations leading to 
multinucleon ejection. The contribution is absent in a free 
nucleon neutrino reaction and in the MB event selection is 
treated as CCQE giving rise to effective large Ma value. 

A microscopic evaluation of the multinucleon ejection 
contribution was reported in [7]. The computations were 
done in the theoretical scheme which was succesfull in de- 
scribing electron scattering in the kincmatical region of 
QE and A peaks together with the dip region between 
them. The model was applied to MB 2D cross section 
data and a fit to the axial mass value was done. In the fit- 
ting procedure [5] the authors included an overall 10.7% 
normalization error. The two-parameter fit gave results: 
Ma = 1.077 ± 0.027 GeV and for the normalization scale: 
A = 0.917 ± 0.029. Using the low-momentum cut proce- 
dure, as proposed in [5], with q cut = 400 MeV the value 
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Contributions to x 2 : Transverse Enhancement Model 



M H 



□ 

'F □ ■ 

□ 
'hp □ 

200 400 



= = □ □ 

a a □ □ □ 
• = □ □ □ 

= □ I I □ a 
■ ■ 

n n ■ 

□ - 
□ 



600 800 1000 1200 1400 1600 1800 2000 
Muon kinetic energy [MeV] 



Contributions to % 2 . : Effective Large Axial Mass Model (M a =1350 MeV) 
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Fig. 1. TEM: contributions do Xmin- Contributions in bins 
are proportional to the area of boxes. If a box is crossed the 
model prediction is larger then the experimental data 



Fig. 3. ELAMM: contributions do \ 
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Contributions to x 2 . : Transverse Enhancement Model (with FG) 
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Fig. 2. TEM-FG: contributions do \l 



M A = 1-007 ± 0.034 GeV was obtained. There is still 
another approach to include 2p-2h contribution coming 
from Meson Exchange Current (MEC) diagrams jTU] . It is 
shown that with the MEC contribution one gets closer to 
the MB experimental results. 



2 Transverse Enhancement Model (TEM) 

In [TT] a new approach to model CCQE scattering on nu- 
clear targets was proposed. The approach is intended to 
be easy to implement in MC event generators. It assumes 
that it is sufficient to describe properly an enhancement 
of the transverse electron QE response function keeping 
all other ingredients as in the free nucleon target case. 

The authors of [TT] proposed a transverse enhance- 
ment function for the carbon target. For low Q 2 its form 
is determined by the scaling arguments while for high Q 2 
(> 0.5 GeV 2 ) it is obtained as a fit to the inclusive electron 
cross section data from the JUPITER experiment. The 



prescription to include TE contribution in the numerical 
computations amounts to rescaling of the magnetic proton 
and neutron form factors: 



G P M n (Q 2 ) -> \/l + AQ 2 eM-^)G p ^(Q 



D 



(2) 



where A = 6 GeV" 2 and B = 0.34 GeV 2 . 

TEM model offers a chance to explain an apparent 
contradiction between recent low (MB) and high (NO- 
MAD) neutrino energy Ma measurements: for energies 
up to ~ 700 MeV the model predicts the CCQE cross sec- 
tion similar to effective large axial mass predictions with 
Ma = 1.3 GeV. For higher neutrino energies the TEM 
cross section becomes smaller and at E v ~ 5 GeV corre- 
sponds to Ma ~ 1.15 GeV. 

The aim of this paper is to confront the predictions of 
the TEM model with the MB CCQE data. 



3 Results and discussion 

In our numerical analysis we compare predictions from 
two models: 

A) effective large axial mass model (ELAMM) wih 
Ma = 1-35 GeV together with the FG model with pa- 
rameter values as in the MB experimental analysis: pp = 
220 MeV and B = 34 MeV. 

B) TEM with the standard axial mass M A = 1-014 GeV 
(as used in |llj). We investigate two implementations of 
the TEM: (i) as in the original paper: without Fermi mo- 
tion and with Pauli blocking effect introduced by means of 
the NEUGEN Q 2 dependent reduction function; (ii) with 
the Fermi motion and Pauli blocking implemented via the 
FG model; we call the model: TEM-FG. 

In both models modifications of the standard (Ma ~ 
1.03 GeV) theory are introduced in the Q 2 dependent way 
in agreement with the MB analysis of the 2D distribution 
of final muons (see Figs 11,12 in [3]). 

We produced three samples of 10 6 events using NuWro 
MC event generator [T^] . We checked that statistical fluc- 
tuations are small. Because in the MB data there is a large 
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Fig. 4. Differential cross section in Q 2 



Table 1. Results for 2D and Q 2 fits 



model 


A2D 


2 

XlD ,min 


A Q 2 




ELAMM 


1.03 


34.1 


1.075 


15.8 


TEM 


1.03 


196.2 


1.015 


22.3 


TEM-FG 


1.135 


133.3 


1.08 


44.0 



overall flux (normalization) error we introduce a renormal- 
ization factor to the x 2 statistical test denned as [8 : 



X 2 (A) 



A" 



AX 



(3) 
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\dT^dc 



V 



A 
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\ dT^d cos 6 J , 



with AX = 0.107. It means that basically we compare the 
shapes of two-dimensional distributions of events. 

Results are shown in Table Q] in the second and third 
columns. The number of degrees of freedom is DOF = 137 
(the number of non-zero bins) -1 = 136. We see that both 

XTEM,min an d Xt E AI — FG ,rnin are la r g er then Xe L AM M ,min ' 



Additionally, the value 



Xelam M /, 
DOF 



34.1 
136 



0.25 is much 



smaller then 1 which suggest that the shape errors were 
evaluated in a too conservative way. 

Figs [TJ [5] and [3] show contributions to Xmin horn three 
models. Contributions in bins are proportional to the area 
of boxes. If a box is crossed the model prediction is larger 
then the experimental data. We see that three patterns 
are rather different in shape and in the case of Figs. [1] and 
[5] they may indicate that the models do not reproduce the 
cross section Q 2 dependence very well. We will come back 
to this point later. In all the cases there is a significant 
deficit of events in the region cos 9^ ~ 1. These are low Q 2 
events for which it is known that techniques going beyond 
the FG (like RPA or CRPA) should be used [13]. 

We made a similar statistical analysis with the Qq E 
differential cross section data. The number of bins is smaller 
(17) but the uncorrelated relative shape errors are also 
much smaller. The \ 2 is defined as: 
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Contributions to x 2 . : Transverse Enhancement Model 
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Fig. 5. TEM: contributions do Xmin with the low momentum 
cut q cu t = 400 MeV/c. 
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Results are shown in Table Q] in the last two columns. 
The number of degrees of freedom is now DOF = 17 
(the number of bins) -1 = 16. The results indicate that 
also here the ELAMM performance is better then that 
of TEM and TEM-FG. Fig g] shows the MB data and 
three models predictions at the best fit points. In all the 
models there is a significant disagreement with the data at 
Qqe < 0- 1 GeV 2 . At large Q 2 QE the ELAMM cross section 
is bigger then that of TEM and TEM-FG in agreement 
with Figs. 7 and 8 from [11] . 

We investigated an impact of the low Q 2 bins on the 
final results and aplied the low momentum transfer cut 
Qcut — 400 MeV/c, as explained in 9 . 2D best fit values 
of x 2 were reduced to 170.8, 98.7 and 23.3 for TEH, TEH- 
FG and ELAMM respectively and corresponding best fit 
values of A were only slightly increased (by about 2-3%). 
Fig. [5] shows the contributions to Xtem min with the low 
momentum cut applied. 

Finally, we tried to understand why for the TEM model 
Xq2 min is relatively small even if there are 2D bins pro- 
ducing very large contribution to \\ D min . We selected 



three such bins: T M e (500, 600) MeV and cos 0^ e (-0.1,0.2) 
with contributions: 9.5, 7.2 and 5.5 (see: [5]). We checked 
that according to TEM events from the selected bins con- 
tribute only to two Q 2 QE bins: (1.2, 1.5) GeV 2 and (1.5, 2) GeV 2 
(cross sections 5.42-10 -48 cm 2 and 1.44-10~ 47 cm 2 /MeV 2 
respectively). Fig [5] shows that these contributions repre- 
sent a small fraction of the overall cross section in the last 
Qqe bin and the large disagreement in three selected bins 
is hidden in the averaged Qqe analysis. It is clear that in 
order to get a deep insight into CCQE, the complete 2D 
data should be analyzed. 

We conlude that it seems that the effective large axial 
mass model leads to better agteement with the MB data. 
However, one should remember that a real challenge is to 
provide predictions for the hadronic final states in mult- 
inucleon ejection effectively described by either ELAMM 
or TEM. 
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